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£n^j ■ Abstract. In the previous paper, we construct the angular distribution functions for muon and electron 

as well as their relative fluctuation functions to find suitable discrimination procedure between muon 
and electron in Superkamiokande experiment. In the present paper, we are able to discriminate muons 
from electrons in Fully Contained Events with a probability of error of less than several %. At the same 
time, our geometrical reconstruction procedure, considering only the ring-like structure of the Cherenkov 
image, gives an unsatisfactory resolution for 1 GeV e and /x, with a mean vertex position error, Sr, 
of 5-10 m and a mean directional error, 50, of about 6°-20°. In contrast, a geometrical reconstruction 
procedure utilizing the full image and using a detailed approximation of the event angular distribution 
works much better: for a 1 GeV e, 5r ~2m and 50 ^3°; for a 1 GeV /x, 5r ~3m and 50 ^5°. At 5 GeV, the 
corresponding values are ~1.4m and ^2° for e and ~2.9m and ^4.3° for \i. The numerical values depend 
on a single PMT contribution threshold. The values quoted above are the minima with respect to this 
threshold. Even the methodologically correct approach we have adopted, based on detailed simulations 
using closer approximations than those adopted in the SK analysis, cannot reproduce the accuracies for 
particle discrimination, momentum resolution, interaction vertex location, and angular resolution obtained 
by the SK simulations, suggesting the assumptions in these may be inadequate. 
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PACS. 13.15.+g Neutrino interactions - 14.60.-z leptons 

1 Introduction First, one needs to calculate patterns for classes \i and e 

with the help of models of Cherenkov light spatial-angular 
In the preceeding paper [lj, we have proposed a new dis- distribution. Second, one should compare the event under 
crimination procedure between electron and muon in the consideration with both patterns by calculating (7-a) 
SK experiment instead of the standard SK procedure, and q el (7.b) which will be shown later (See page 2). Third, 
For the purpose, we have constructed the mean angu- one compares their difference with q crit thus making a 
lar distribution functions for the charged particles con- decision on the type of event. 

cerned (muons and electrons/positrons), as well as the In order to discriminate electron neutrino events from 
corresponding functions for the relative fluctuation. In the muon neutrino events, we construct discrimination func- 
present paper, we apply our discrimination procedure to tions in the following way. 

the SK experiment. We consider each event to be a random vector Q = 

{Qj} containing the contributions to all photomultipliers 
(PMTs) of the detector. Here j is the PMT index: j = 

2 Procedure for discriminating electron ^' ""' ^ anc * ^ * s tot ^ num ber of PMT of the detec- 

° tor. Inere are two classes ol events to be considered by the 

neutrinos from muon neutrinos: Event type type identification procedure: uo x = e (electron initiated) 

definition and U2 = /i (muon initiated). Monte Carlo simulations 

of event optical images provide a possibility to study the 

Figure 1 shows the process of classification, i.e. how event behaviour of images of both classes, namely, the image dis- 

is called mu-like (muon event) or e-like (electron event), tribution functions F{Q X , Q 2 , Qn] ^i, E , r , o ) which 

are simultaneous distribution functions of PMT contribu- 

Send offprint requests to: tions Qj for given particle type Cc^, energy Eq, injection 
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Fig. 1. schematic view of pattern recognition. 

point ro and direction Oo. i.e. the most differential and 
complete characteristics of the classes. 

It is, however, too expensive to deal with these func- 
tions: in principle, one should simulate an enormous sam- 
ple to know the most differential distribution function 
even for only one set of primary parameters {^, Eq, r , Oo}. 
A more practical solution is to determine only the mean 
vector {Qj((Ji, Eq, ro, Oo)} and the covariance matrix 
UQ(uJi, Eq, ro, 0o) = cov(Qj,Q m ) and then assume the 
distribution is a multi-variate normal one: 



p(Q;Wi,Eo,r ,e ) = (2w)- N / 2 ■ (detS- 1 
xexpl-iQ-Q^fs- 1 (Q-Q"" 



-1/2 



(i) 



Here 

Q W1 = Q(^,^ ,r ,6>o),^ = ^Q(^,^o,r ,6> ). 

In this case, one still needs a sufficiently large sample to 

calculate Q and Eq accurately. 

Our approach assumes that one can calculate the mean 
image vector Q and fluctuation vector SQ using the ap- 
proximations of the Cherenkov light mean angular distri- 
butions and their fluctuations. Simulation shows that cor- 
relation coefficients of single PMT contributions amount 
to 0.6(0.8) for the neighbouring PMTs and go down to 
0.1(0.1) for distant PMTs in electron(muon) events, which 
makes the covariance matrix diagonal-like. 

Thus, we neglect all correlations between Qj and the 
resulting normal distribution becomes simpler: 



p(Q;cj;,£;o,ro,#o) = (2tt) 



-N/2 m 




xexp • 



i=i 



SQT 



, (2) 



For a certain event geometry (starting point r and 
direction 0q of the particle) and energy Eq we calculate 



mean pattern images = (Eq^yq^Oq) and their 
deviations 5Q^ Z = 5Q^(Eq, r , o ) : 

Q^(E ,ro,0o) 



^2 — • C0S Xj,k • exp 

k=l P 3* 



^abs 



(3) 



k=l 



s 

— • cosxj.k • exp 

P J,k 



Pj,k 
^abs 



xK'"(0i,k)-<sn*i,k)] . 



(4) 



where k is the segment index, n is the number of segments 
in a track/shower, S is the PMT area, pj : k is the distance 
between a segment center and a PMT center (see Fig. 2), 
Xj,k is the cosine of the angle between pj^ and a PMT axis 
(which is normal to the tank surface), 0j^ is the radiation 
angle from the segment center to the PMT center, and 
Xabs is light absorption length in water. And F^{0j^) 
is the mean angular distribution function, S^(0j^) is its 
relative fluctuation. 

Eq.(4) implies there are no correlations between the 
contribution from different segments of track/shower which 
is close to reality: in simulated electron event samples cor- 
relation coefficients between neighbouring segments do not 
exceed 0.4 and amount to 0.1 for distant segments. In 
muon events characteristic values of the coefficients be- 
tween different segments (including neighbouring ones) 
are even smaller: 0.1; the exceptions are the trailing (last) 
segments which emit about 100 times less than the other 
ones but are tightly correlated due to the muon decay 
process. 

Then we apply Bayes' decision rule which minimizes 
the error of misidentification assuming that prior proba- 
bilities of both types of events are equal: 

P(e/Q) _ p(Q/e) 

r ~ P(/VQ) ~ p(Q/a0 ~ 



(5) 




N 



1/2 



exp < 




-Qj) 2 m} 


exp | 


N 

- E (Qj - 

k 3=1 


-Qj) 2 / s Qj^j 



where p (Q/e) and p (Q//i) come from Eq. (2) . 

The simplest criterion, q, which we use to define the 
type of event is: 



q = Inr = - q ei + C , 



(6) 



N 



= T,&r -q?) / 5 q?> ( 7 -^) 

3 = 1 
N 

Qei = £(Qf p ~ Qjf/ 5 Qj > ( 7 - & ) 
j=i 
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Fig. 2. geometry for eqs.3 and 4. 
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i J}i SQBj 

3=1 



(7- 



where Qj Xp is the contribution to j-th PMT in the "ex- 
perimental" image under consideration. The event is con- 
sidered to be e-like if q > and /i-like if q < 0. 

A slightly more general form of the criterion can help 
to reduce the errors of type definition: 

q = q a - A ■ q el + B , (8) 

where A, B are tuned to minimize the identification errors; 
the usage is the same as above. 



3 Examination of the procedure for 
discriminating electrons from muons 

3.1 A virtual 1 kilo-ton water tank detector for the 
KEK experiment 

As explained in the preceeding paper [I], the SK discrim- 
ination procedure was constructed on the following as- 
sumptions in the Kamiokande analysis - the predecessor 
to SK : 



N^ exp (direct) = a e x N M c(0i,Pe) x 



xcxp |~|x /(©), (9) 



16.9 



N iyexp (direct) = 



a,, x 



Utsintii + U x (£)) 



x sin 2 0i + N iiknock (0i) > x 



xcxp | j x /(©), (10) 



p / 1 ( ( N obs-N exp y 

Pr oo = ■ exp ' 



2trt 



2a 2 



Ppattern (^) 



exp 



pattern 



1 f X 2 (e) - min[x 2 (e),x 2 (l^)r 2 



(11) 



, (12 -a) 



(M) 



exp • 



X 2 (/i) -min[x 2 {e),x 2 {lA] \ 
cr x 2 J 



(12 - b) 



Pangle{&) — 

= constant x exp • 

Pangle{f^) = 

= constant x exp < 




(13 -a) 



'obs 



AO 



(13-6) 



^P(e) = Ppattern(e) x P ang i e (e) , (14 - a) 

P(/i) = Ppattern(v) X P a ngle(^)' (14 - 6) 

These expressions was given firstly in Takita [2]. 

However, in order to show the validity of this discrim- 
ination procedure which had never been confirmed by ac- 
cerelator beam, the SK group performed an accelerator 
experiment using the KEK 12 GeV proton synchrotron 
(Kasuga[3j, Kasuga et al.[4j, Sakai[5j, Kasuga[6j). For the 
confirmation, one kiloton water Cherenkov detector was 
constructed and the SK group studied the particle identi- 
fication capabilities over the momentum range 250 MeV/c 
to 1000 MeV/c for muons and 100 MeV/c to 1000 MeV for 
electrons. They concluded that the particle misidentifica- 
tion probabilities with SK discrimination procedure is less 
than 1 to 5% (Sakai, p. 141 [5j, Kasuga, p. 77 [6]). Based on 
these procedures, the SK group have reported experimen- 
tal results which have lead to the claimed existence of neu- 
trino oscillations. To study the positional dependence of 
the discrimination, data was taken for 10 injection points 
in the 1-kiloton detector: WUS, ESS, WUC, WUD, ... (See 
Sakai [5] in detail). We have constructed a virtual 1 kilo- 
ton water Cherenkov detector, including the same PMT 
configuration, as the SK test one in our computer. In the 
virtual detector we simulate physical events and examine 
the influence of the virtual physical events on the virtual 
detector. 

In Figure 3 we give the relation between the emit- 
ted Cherenkov light and the detected Cherenkov light for 
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Fig. 3. The average behavior of the Cherenkov light for both 
emitted and detected as the functions of momentum for muons 
and electrons which are generated at WDS. 



the WDS injection point (see Sakai's paper). The emitted 
Cherenkov light denotes the total Cherenkov light pro- 
duced by the charged particle concerned, while the de- 
tected light records only that which is measured by the 
detector. This relation holds almost exactly, irrespective 
of the generation point, because the absorption length is 
so large compared with the dimensions of the 1 kilo-ton 
detector. 

It is clear from Figure 3 that there is big difference 
as for the Cherenkov light quantity between the emit- 
ted Cherenkov light and the detected one. The emitted 
Cherenkov light denotes, by the definition , the total Che- 
renkov light, while the detected Cherenkov light denotes 
only small part of the emitted one. Consequently, we ex- 
pect variety of the patterns of the detected Cherenkov 
lights for the events concerned due to the PMT configu- 
ration. 

Using the discrimination procedure developed in sec- 
tion 2 (Eq.(8), A=l, B=q cr it), we try to discriminate 
300 MeV electrons from 500 MeV muons, both of which 
result, on average, in the production of the same amount 
of Cherenkov light. The result is given in Figure 4. As can 
be seen, discriminating between the electron and muon 
seems to be very difficult. In Table 1 we give misidentifi- 
cation probability for different q cr it- The parameter q cr u 
denotes a criterion for discriminating between muons and 
electrons: e — > \i denotes the probability for an electron 
to be misidentified as a muon, and \i — > e denotes the 
probability for a muon to be misidentified as an electron. 
It is easily understood from Figure 4 that a non-negligible 
number of muon events penetrate into the electron events 



Table 1. Misidentification probability between 300MeV elec- 
tron events and 500MeV muon events in the virtual KEK 1 
kiloton tank detector, taking into account forward scattering 
only. 



run 


Qcrit 






1 


-4.00e2 


1.67 


25.67 


2 


-5.00e2 


2.67 


22.33 


3 


-6.00e2 


6.67 


19.67 


4 


-7.00e2 


12.67 


18.67 


5 


-8.00e2 


20.33 


17.30 


6 


-9.10e2 


27.33 


15.67 


7 


-1.00e3 


35.67 


12.67 



KEK tank, WUS, 300*0.3GeV e + 300*0.5 \i 
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Fig. 4. The q e i-q^ distribution for simulated events for 300 
MeV electrons and 500 Mev muons in the KEK kiloton tank 
with angular distribution function describing only the forward 
hemisphere. These particles produce approximately the same 
amount of Cherenkov light. The notation " 300*0.3 GeV e " 
denotes the reslutls for 300 MeV electrons from 300 simula- 
tions. 



region, so that some muon events may be misidentified as 
electrons. 



Thus, we carefully examined the muon events which 
were located in the electron event parameter space, i.e., 
where electron events are found with high probabilities. As 
a result, we found that the majority of such muon events 
are fi — e decay events which should be identified as muon 
events. However, from the view point of the pattern recog- 
nition, we could not recognize whether the events in the 
gray zone between the 'clearly electron events' and the 
'clearly muon events' were electrons or muons. To resolve 
this, we extend our angular distribution function, which 
was defined only within the forward hemisphere, to in- 
clude the backward hemisphere and redefine q e \ and q^ in 
accordance with this extension. 
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Table 2. Misidentification probability between 300MeV elec- 
tron events and 500MeV muon events in the virtual KEK 1 
kiloton tank detector, taking into account forward and back- 
ward scattering. 



run 


Qcrit 






1 


2.00e2 


0.00 


4.67 


2 


1.00e2 


0.33 


3.67 


3 


0.00 


0.67 


3.00 


4 


-1.00e2 


0.67 


2.00 


5 


-2.00e2 


2.00 


2.00 


6 


-3.10e2 


3.67 


2.00 


7 


-4.00e2 


12.00 


2.00 


8 


-5.00e2 


20.33 


1.67 


9 


-6.00e2 


30.00 


1.67 


10 


-7.00e2 


43.67 


1.33 



KEK tank, WUS, 300*0.3GeV e + 300*0.5 \i 
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Fig. 5. The q e i-q^ distribution for simulated events for 300 
MeV electrons and 500 Mev muons in the KEK kiloton tank 
with angular distribution function describing both forward and 
backward hemisphere. These particles produce approximately 
the same amount of Cherenkov light. The notations are the 
same as for Figure 4. 



After such a manipulation, one gets much smaller con- 
tamination of electron domain by muon events. 

In Figure 5, we give the revised q e i-q^ distribution 
where there is almost no intersection between the muon 
events and electron events. In Table 2, we give the misiden- 
tification probability based on the revised procedure. Com- 
paring Table 2 with Table 1, we find that misidentification 
probability has decreased remarkably. 



The probability distribution of the standard SK esti- 
mator for particle identification, if it were assumed to be 
remotely comparable, should be compared to Figure 4, not 
Figure 5, because the effect of fi — e decay was not consid- 
ered correctly in the SK simulation (See Sakai,p.l54, [5]). 



There are, thus, big differences between our result and the 
SK results (For example see Figure 3 in Sakai[5j). Namely, 
our initial result (Figure 4) could not reliably separate 
electron events from muon events, although it is claimed 
the standard SK analysis separates electron events from 
muon events almost completely. We can not imagine such 
a complete discrimination could be realized by the SK 
procedure. 

The big difference between the "complete" discrimi- 
nation by the SK (Figure 3 in Kasuga et al.[4j) and the 
poorer discrimination obtained by us (Figure 4) extends 
further, if we take into consideration the uncertainties 
added by the real detector to our Figure 4. This is because 
we do not consider the effect of photoelectrons produced 
from the parent Cherenkov light, which obeys a stochas- 
tic physical process, or the effect of scattered Cherenkov 
light. These two factors must increase the uncertainties in 
the discrimination between muons and electrons. 



3.2 Discrimination of electrons from muons in the 
Super-Kamiokande detector in our simulation 

Also, we have constructed a virtual SK detector incorpo- 
rating the actual PMT configuration and have carried out 
the experiment numerically in our virtual SK detector. 
Based on the revised procedure for the discrimination of 
electron events from muon events developed in the pre- 
vious section, we have discriminated electron events from 
muon events. 

In Figure 6 we give the discrimination of the muon 
events from electron events in the virtual SK detector by 
the revised method. We conclude that we could discrim- 
inate muon events from electron events using the revised 
method with high probability. In Figure 7 we show the 
scatter plot between q e and q^. From Figure 6 and Fig- 
ure 7, we can conclude that our procedure for the discrim- 
ination of electron events from muon events works well. 

In Figure 8 and Figure 10 we show the discrimination 
of muon events from electron events for a comparison of 
500 MeV electrons with 700 MeV muons, and 700 MeV 
electrons with 900 MeV muons, respectively. These elec- 
tron and muon energies were chosen as they yield roughly 
the same amount of Cherenkov light. 

Also, scatter plots corresponding to Figure 8 and Fig- 
ure 10 are given in Figure 9 and Figure 11, respectively. 
Comparing Figure 6 with Figure 8, it is apparent that the 
latter allows a better discrimination. Comparing Figure 8 
with Figure 10, the latter is also better than the former 
for discrimination. The reason for this is that the gener- 
ation of Cherenkov light for muon events becomes more 
different from that of the electron as their primary ener- 
gies increase. Discrimination procedures of muon events 
from electron events in the SK detector are examined for 
events whose primary particles start at the "scaled WUS" 
position (#0=986 cm, yo=0 cm, zo=282 cm, = 7r/2and 

= o). 
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Fig. 6. The q e i-q^i distribution for simulated events for 300 
MeV electrons and 500 MeV muons in the virtual SK detector. 
The notations are the same as for Figure 4. 



Fig. 8. The q & i-q^ distribution for simulated events for 500 
MeV electrons and 700 MeV muons in the virtual SK detector. 
The notations are the same as for Figure 4. 
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Fig. 7. The scatter plot for q e i and q^ for 300 MeV electrons Fig. 9. The scatter plot for q e \ and q^ for 500 MeV electrons 
and 500 MeV muons considered in Figure 6. and 700 MeV muons considered in Figure 8. 
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Fig. 10. The q e \-q^ distribution for simulated events for 700 
MeV electrons and 900 MeV muons in the virtual SK detector. 
The notations are the same as for Figure 4. 
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Fig. 11. The scatter plot for q e \ and for 700 MeV electrons 
and 900 MeV muons considered in Figure 10. 



4 Error Distribution for the vertex point and 
direction for muon and electron events 

4.1 TDC procedure 

The SK analysis introduces a "TDC procedure" to deter- 
mine the vertex position. The principle of the TDC proce- 
dure is to find the position where the time residuals, for 



the PMTs being fitted are minimized. The time residual 
tj of the i-th PMT is denned as 



■xy/(x-Xi)* + (y- yi ) 2 + (z-z^ , (15) 



where t® is the hit time of the i-th PMT, (xi,yi,Zi) is 
the position of the i-th PMT, (x, y, z) is the effective 
emitting point position and c/n is the velocity of the 
Cherenkov light in water. That all the light is emitted 
from the same effective point in space and comes to j-th 
PMT exactly at the moment corresponding to the mean 
time tj of the j-th PMT Cherenkov pulse is not really 
true. A simple equation for time residual used in a x 2 ~ 
like sum (the system of linear equations for the effective 
point coordinates (x*,y*,z*) is overdefined!) can give ef- 
fective point estimate after the sum's minimization (with 
respect to (x*,y*,z*)) even if the original assumption is 
not valid. The effective point thus deduced does not co- 
incide with the center-of mass of the light emitting sys- 
tem (e-shower or //-track) because of specific mechanism 
of Cherenkov pulse formation and will usually differ from 
the event starting (injection) point. 

The TDC procedure based only on tj cannot estimate 
the event direction because to define a direction one needs 
at least two points. Direction estimates could be obtained 
as a result of Cherenkov pulse shape analysis for each suf- 
ficiently illuminated optical module if the PMT and elec- 
tronics are fast enough for such analysis. 

Sakai shows the time residual distribution of typical 
event (a 1 GeV/c, electron) which is distributed over 50 
nanoseconds (Sakai, p. 38 [5j), assuming a point-like source. 
We simulate the Cherenkov light in the cascade shower 
using GEANT 3.21 and the tools we have developed. In 
Figure 12 we give one example for the time residual distri- 
bution for a 1 GeV primary electron based on Eq.(15) with 
the use of the detailed simulation of the cascade shower. 
In our calculation, we simulate shower particles and the 
accompanying Cherenkov light due to shower particle con- 
cerned in an exact way. Then, we know the starting point 
of the primary electron. Shifting the starting point from 
the real point to a range of artificial ones, we can obtain 
the time residual distribution for each position, and ex- 
amples are given in Figure 12. 

Among the five different starting points, which includes 
the true one, the smallest standard deviation is obtained 
in the case of Z' — Z + 100 (cm), where Z' and Z denote 
the assumed vertex point and the real vertex point, respec- 
tively. Thus, the apparently most probable vertex point is 
not real one, but is offset by 100 cm from the real one. Of 
course, this is only one example and not average behavior. 
However we examined many individual cases and confirm 
that this is usual character which should hold even for the 
average behavior. 

The comparison of our simulations with the experi- 
mental data from the SK experiment (Sakai) reveals a 
large difference. The width of our time residuals distribu- 
tion is within one nanosecond, while the width for SK is 
~50 nanoseconds. The probable reasons are that we have 
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Fig. 12. One example of the time residual distributions for 1 GeV electron primary shower, assuming different vertex positions. 



not considered the PMT and electronics response func- 
tions, and that we have neglected light scattering. 



We calculate the time residuals for electrons of 500 MeV, 
1 GeV, 3 GeV and 5 GeV, assuming that all Cherenkov 
light comes from certain points of a shower/track. From 
these calculations one can see that the smallest time resid- 
uals do not give the vertex position but yield points shifted 
from the vertex point along the direction of the cascade 
shower, namely, 50 to 100 cm for 500 MeV electrons, 100 
to 150 cm for 1 GeV (FigurdT2|) and 3 GeV electrons, and 
150 to 200 cm for 5 GeV electrons. Such a tendency is quite 
understandable if we consider the size of the shower/track: 
the effective point should not be too close to the starting 
or ending points. The error of effective point location by 
minimizing the width of the distribution amounts to about 
50 cm. 



From the much larger width of the SK time residual 
distribution it is clear that in experimental conditions the 
effective point location error should be a few times greater 
because the minimum of the width as a function of effec- 
tive point position would be much less pronounced. 



For reasons mentioned above, we conclude that the 
SK TDC procedure is not suitable for the determination 
of an accurate vertex position for electron events. The sit- 
uation for the muon events is essentially the same but 
must be worse than in the case of electron events as muon 
events have a longer extent than the corresponding elec- 
tron events. 



4.2 ADC procedure 

The standard SK analysis uses an "ADC procedure" for 
event geometry reconstruction which is based only on the 
total amount of Cherenkov light detected by the PMTs. 

We have constructed and used an ADC procedure which 
is based on the more detailed model for the angular dis- 
tribution for the Cherenkov light described above. 



4.2.1 Fundamental procedure for the error distribution for 
muon events and electron events: event geometry definition 
procedure 

Our geometry reconstruction procedure assumes the type 
and energy of event are known. Thus, we only need to 
determine the event starting point r and the direction of 
the primary particle 0q. 

The first step is to find a ring-like (or arc-like) struc- 
ture in the image. In muon events we observe two ring- 
like structures because the muon-decay electron produces 
a ring-like or a spot-like structure of its own, however, we 
do not consider muon energies below 500 MeV, so in muon 
events we choose the most intensive ring-like structure. 
This task is accomplished by row-by-row (or column-by- 
column) scanning of the arrays, representing an image, for 
maxima above a certain threshold. 

The second step is defining the first approximation 
i*i , #i of geometrical parameters by fitting a ring struc- 
ture with the following simple (cone-like) model: all light 
Qtot emitted by a /i-track/e-shower is assumed to come 
from one effective point W on the track with the angular 
distribution of light being Heaviside-like: 
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cone-like light angular distribution 




Fig. 13. 'cone' figure. 



detailed light spatial-angular distrit 




F(0) 



Fig. 14. 'full' figure. 
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sinOdO = Q tot (16) 



The zero approximation for geometry parameters could 
be arbitrary because the ring-like structure is usually very 
pronounced and we use the cone-like model of Cherenkov 
light angular distribution with abrupt edges, the function 
to minimize with respect to r, is: 



Al 



G(r,0) = J2 



(QrvoA) 



1=1 



Qi xp (ro,0o) 



(17) 



where I is the PMT index within the ring structure, M is 
the number of PMTs in the ring structure, Q z ercp (ro, #o) is 
the contribution to Ith PMT from the event under consid- 
eration, Qi(r,9) is an estimate for this contribution cal- 
culated with the cone-like angular distribution model de- 
scribed above, F{0): 



Qi(r,0) = 



S 
Pi,w 



• cosxi,w ' exp 



Piyv 

^abs 



■F(8i,w), (18) 



The third and final step is the improvement of geome- 
try parameter estimates by fitting all of the image above 
a certain threshold, Qthr, with the pattern image 
{Qj(Eo, r, 0), 5Qj(Eo, r, 0)} calculated using the detailed 
model of Cherenkov light angular distribution. The previ- 
ous step fit result ri , 1 is used as the first approxima- 
tion, and the function to minimize with respect to r, is: 



H{v,e) 



r-Q- xp >Qth, 



(Q; xp (E o ,t o ,0 o ) 



Qj{E ,v,e)Y 



5Qj(E ,r,e) 



(19) 



The PMT contribution threshold, Qthr, could be op- 
timized with respect to the resulting geometry parame- 
ter resolution, as we show later. The optimal threshold, 
Qthr,opt, increases with primary energy, which could be 
used to enhance the resolution, because an estimate of the 
primary energy can be obtained from the total amount of 
Cherenkov light detected quite independently of the ADC 
procedure. 

Figure 13 shows the process of geometry reconstruc- 
tion using ring structure. Assuming geometry parameters 
of an event to be r, 0, cone-like angular distribution of 
light (Eq.(16)) is used to evaluate contributions (18) to 
all detectors (PMTs) of the ring structure in order to cal- 
culate G(r, 0) (Eq.(17)) which compares this rough light 
pattern with the light ring structure. Adjusting r and 
(see arrows in the figure) so that G approaches its mini- 
mum one obtains a rough estimate of the true geometry 
parameters ro, 0o. 

Figure 14 shows the process of geometry reconstruc- 
tion using full light image of an event. Assuming geometry 
parameters of an event to be r, (for instance, these are 
values obtained in ring structure processing), the model 
of spatial-angular distribution of light is used to evaluate 
pattern image {Qj{E^ 1 r, #), 5Qj(Eo, r, 0)} which is then 
used to calculate i?(r, 0) (Eq.(19)) which compares it with 
the full light image of the event. Adjusting r and (see ar- 
rows in the figure) so that H approaches its minimum one 
obtains a better estimate of the true geometry parameters 



4.2.2 Application of the technique to simulated events 

By using the technique developed in the previous section, 
we analyze simulated events to determine the error distri- 
butions for the vertex position and for the particle direc- 
tion. Here, we examine the error distributions for 300 MeV 
electrons and 500 MeV muons, which yield roughly the 
same amount of Cherenkov light. 

In Figure 15, we give the error distribution for the ver- 
tex position for 300 MeV electrons for different Cherenkov 
threshold quantities. "Ring only" denotes that only the 
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information from PMTs whose Cherenkov photons con- 
tribute to the Cherenkov ring are used for the estimation 
of the error. "Full proc., thr=lph" denotes that informa- 
tion from "ring only" PMTs and also those exceeding 1 
Cherenkov photon are utilized for the estimation on the 
error. For "full proc, thr=5" and "full proc, thr=10", this 
latter threshold is raised to 5 and 10 photons, respectively. 

In Figure 16 the error distribution for the vertex posi- 
tion for 500 MeV muons are given. It is clear that a wider 
error distribution is obtained for a "ring only" analysis. 
A narrower error distribution results from the "full proc, 
thr=10 ph" algorithm. This is the same as in the case of 
the electron. However, muons generally have wider error 
distributions than electrons: the mean error for 500 MeV 
muons for the vertex determination in the full analysis is 
2.9 m while for 300 MeV electrons it amounts to 2 m. 

In Figure 17, the error distribution for the direction of 
the 300 MeV electron is given. As expected, "ring only" 
gives the largest uncertainty distribution, while "full proc, 
thr=10ph" has the narrowest error distribution, with a 
mean error of about 3.7°. In Figure 18, we give the corre- 
sponding distributions for muons. The same trend is seen 
as for electrons, though the muons have a wider uncer- 



tainty distribution. The mean direction uncertainty in the 
best case is 4.9°. 

Now, we compare 1 GeV electrons with 1 GeV muons, 
both of which yield roughly the same quantity of Cherenkov 
light. In Figure 19, we give the error distribution for the 
vertex position for 1 GeV electrons in the case of "full 
proc,thr=10" . The mean error is 1.9 m. In Figure 20, the 
corresponding quantities for the muon are plotted. The 
average error for the vertex position is 3.2 m. Again, the 
error of the vertex point for muons is larger than for elec- 
trons. 

In Figure 21, the error distribution for the direction 
for 1 GeV electron is shown. The average direction error is 
3.0° for "full proc,thr=10". The corresponding quantities 
for the muon are plotted in Figure 22, where the mean 
error is 5.3°. Once again, the directional error for muons 
is larger than that for electrons. 



In Table 3, we summarize the error distributions for 
the vertex points and the direction for both electrons and 
muons. Both mean errors and root mean square errors 
are given. Errors are also given for the different criteria, 
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Table 3. Mean and standard deviation of the error in the vertex position and the direction due to primary electrons and 
primary muons. These are given for different criteria for the Cherenkov threshold. Ring proc. denotes errors estimated using 
the Cherenkov ring only. [1], [5], [10], [20] denote errors estimated by the combination of [ring proc] with a Cherenkov photon 
threshold of 1, 5, 10 and 20photons, respectively. Aim denotes the mean direction error in degrees. Als denotes the standard 
deviation for the corresponding mean values. Rm denotes the mean position error in metres. Rs denotes the standard deviation 
for the corresponding mean values. 
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Error distribution for the vertex position for 1 GeV Fig. 20. Error distribution for the vertex position 1 GeV 



from the injected scaled WUS. The number of simu- 
100. 



muons from the injected scaled WUS. The number of simula- 
tion is 100. 



namely, different Cherenkov light threshold. From Fig- 
ures 15 to 18 and Figures 19 to 22 and Table 3, it should 
be noticed the followings: 

1. Of the different criteria considered, the "ring only" 
procedure results in the largest error. The reasons are 
as follows: The concept of the ring structure is essen- 



tially fuzzy, both in our procedure and the SK proce- 
dure, and information from ring structure is only part 
of the total information available for the pattern recog- 
nition. It is, therefore, natural that the vertex position 
and directional errors are largest in the "ring only" 
analysis. The standard SK analysis uses ring structure 
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Fig. 22. Error distribution for the direction 1 GeV muons 
from the injected scaled WUS. The number of simulation is 
100. 



only, and their errors are amplified by that fact that 
the analysis ignores fluctuation effects. 

2. Muons events have larger uncertainties than electron 
events. For both electrons and muons, the sources for 
the Cherenkov light are not point-like and have some 
extent in both cases. Significant errors come from the 
point-like approximation for electron events. 

3. The optimal Cherenkov threshold for the third step 
of geometry reconstruction procedure depends on the 
primary energy of the particle concerned. For energies 
less than 1 GeV, third step with 10 photon threshold 
gives the best results for both muons and electrons 
among the alternatives considered. For 5 GeV electrons 
and muons, 20 photon threshold seems to be optimal 
for the third step. 

4. The fact that the uncertainties for the determination of 
the vertex point and direction are rather large comes 
from the effect of fluctuations, namely the nature of 
the stochastic process concerned (an electron cascade 



account correctly. Even if additional errors exist, they 
should be negligible compared to the uncertainty caused 
by fluctuations. The rather large errors for the vertex 
point and the direction obtained by our model could 
not be reduced substantially, reflecting the essential 
nature of the physical processes concerned. 

5. The SK analyses, according to all published accounts, 
completely neglect fluctuations and also use point-like 
approximations for the electron cascade. Moreover, they 
neglect the scattering effects on muon track geometry. 
As we showed earlier, their simple approximations dis- 
tort the mean values in certain parameter domains. 
The most probable reason for their low error estimates 
is the fact that they completely neglect fluctuations in 
the event development. Our results contradict clearly 
the fine positional resolution of 23 to 56 cm claimed by 
for the SK analysis (Kibayashi, p.73[7j). 

6. It should be noticed that errors derived by us are lower 
limits. As already mentioned, we do not consider the 
production of photoelectrons in PMTs, and only con- 
sider direct Cherenkov photons in our discrimination 
procedures neglecting the diffusion of Cherenkov pho- 
tons. If we include these factors in our procedure, then, 
the actual errors for the vertex position and the direc- 
tion should be larger than that given here. 

7. Here, we make remarks to the results obtained by Mit- 
sui et al([8j. They have calculated photoelectrons from 
the Cherenkov light due to both muon-like events and 
electron like events in the SK experiment by full Monte 
Carlo Method and have showed the clear deviation of 
the photoelectrons concerned from the Poisson distri- 
bution which SK assume(See,Eq.(ll) in the present 
paper). Further, they have calculated PID parameter 
values following the SK PID methods which consists of 
the pattern and opening angles of the Cherenkov light. 
As the resluts of it, they have concluded that the total 
misidentification for muon and electron events is larger 
than or equal to 20 percent in sub-Gev and also at sev- 
eral percent in multi-GeV region. They attribute the 
deviation from the Poisson distribtuion to the neglect 
of the fluctuation from the stochastic proceeses inher- 
ent in electron-like eventsand muon-like events. We 
agree to their interpretation. However, we are forced to 
criticize the adoption of the same PID ESTIMATOR 
as the SK's. We have already pointed out the inade- 
quency of SK PID ESTIMATOR ( See Galkin [lj, page 
6). 



5 Summary 

(1) 



The SK TDC procedure 

The TDC procedure assumes that the Cherenkov light 
originates from a point, and thus does not determine the 
vertex position accurately, because the sources for the 
shower or sequence of muon interactions with the mediumpherenkov light have a no n- negligible extent. In order to 
The utility of model developed in this paper, the mov- utilize the TDC meaningfully, we should take into account 
ing point approximation model, is guaranteed, because the extent of the source for the Cherenkov light in space 
it gives mean values and relative fluctuations precisely and time. Further, ideally we should utilize not only ar- 
and takes all necessary geometrical considerations into rival time of the Cherenkov light but also shape of the 
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pulse in the PMTs. 

(2) The SK ADC procedure: the estimator for particle 
identification 

It is no exaggeration to state that the credibility of the 
discrimination of muon events from electron events in the 
SK analysis is entirely dependent on the validity of the 
estimator adopted by the SK. 

The estimator for particle identification was introduced 
into the analysis of the neutrino events in the Kamiokande 
detector (Takita[2j). The main purpose of the KEK 1 kilo- 
ton experiment was to test, using particles of known type 
and momentum, the validity of the Monte Carlo simula- 
tion used in the Kamiokande detector (Kasuga, p.22[3j). 
The results of this test established the validity of the esti- 
mator and thus it was applied to the standard SK analysis 
(Kasuga et al.[4|, Sakai[5], Kasuga [3]). As a result of us- 
ing this estimator, it has been asserted that evidence for 
neutrino oscillations between muons and tauons has been 
obtained. 

As we have pointed out, the SK estimator for particle 
identification consists of three parts. The first is the proba- 
bility function, in which the spatial part is separated from 
the angular part. The second is the probability function 
for photoelectron production. The third is the introduc- 
tion of the point-like approximation into the Cherenkov 
light for electrons. As we have made clear, the first two 
lack a firm theoretical foundation. The third is an over- 
simplification which leads generally to a wrong estimation 
for the quantity of Cherenkov light. Consequently, the SK 
estimator for particle identification can not be reliable for 
an accurate discrimination between electrons and muons. 
SK results depending on the discrimination of electron- 
like events from muon-like events (see Kasuga et al.[4]) 
must be re-evaluated in this light. 

(3) Errors for the vertex point and the direction 

Only the correct estimator for particle identification 
can estimate the error for the vertex point and direction 
quantitatively. Our estimation (Figures 15 to 18 and Fig- 
ures 19 to 22 and Table 3) shows non- negligible and in- 
evitable errors for the vertex position and direction. It 
should be, particularly, noticed that the fluctuations in 
error in both vertex position and direction are too big. 

Kibayashi (p.73[7j) concludes that the uncertainty for 
vertex point is from 23 cm to 56 cm and the uncertainty for 
the direction from 0.9° to 3.0° using both the estimator for 
particle identification and the TDC adopted by the SK. 
As we have demonstrated, these appear to severely un- 
derestimate the error distributions, which is too far from 
reality. 

(4) In the present paper, we do not take into consideration 
photoelectrons produced by the Cherenkov light in the 
PMT, and we neglect scattering of the Cherenkov light. 
Therefore, our results on discrimination of electron events 
from muon events, only yield lower limits to the realisti- 
cally achievable experimental erroes. If we consider these 
factors, the situation becomes much more complicated. 
In Figure 23, we give the Cherenkov light distributions 
for given numbers of photoelectron. It is easily under- 
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number of simulation is 100. 



stood from the figure that the patterns of electron events 
and muon events are much more vague, because there is 
non- negligible distribution of the photoelectron for a given 
number of Cherenkov photons. Further, we can say that 
it is quite sufficient for us to limit our analysis to the level 
of direct Cherenkov light, neglecting scattered Cherenkov 
light, to clarify the essential points in the discrimination 
procedure between electrons and muons. 
(5) As a result of the points mentioned in (4), our results 
give lower limits for the primary parameter uncertainties 
in the problems concerned. We show a clear separation 
between electron and muon by using a procedure based 
on detailed Cherenkov light angular distribution approx- 
imations. Even though we show a clear separation be- 
tween electron and muon for the same quantity of the 
Cherenkov light, this remains at the theoretical level and 
does not provide the optimal experimental means of sepa- 
rating electrons from muons. This problem is beyond the 
scope of this paper and is better undertaken by those with 
a detailed understanding of the specifics of the detector. 

A part of the preceding paper pQ and present paper are 
found in [9]. 
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